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Branching within branching II: Limit theorems 


Gerold Alsmeyer and Soren Grottrup 


Abstract This continues work started in ID on a general branching-within-branching 
model for host-parasite co-evolution. Here we focus on asymptotic results for rele¬ 
vant processes in the case when parasites survive. In particular, limit theorems for 
the processes of contaminated cells and of parasites are established by using martin¬ 
gale theory and the technique of size-biasing. The results for both processes are of 
Kesten-Stigum type by including equivalent integrability conditions for the martin¬ 
gale limits to be positive with positive probability. The case when these conditions 
fail is also studied. For the process of contaminated cells, we show that a proper 
Heyde-Seneta norming exists such that the limit is nondegenerate. 
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1 Introduction 

We start with a brief review of a general branching-within-branching model for 
the evolution of a population of cells containing proliferating parasites. The model 
has been introduced in the companion paper HI, and we refer to this paper for a 
more detailed introduction of the basic branching within branching process (BwBP) 
(see (|^ below), its relation to other branching models and an account of relevant 
literature. 
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Given the infinite Ulam-Harris tree V with root 0, let the cell tree be formed by 
the subfamily (A^v)vgV of independent Np-valued random variables with common 
law {pk)k>o and finite mean v, which is a standard Galton-Watson tree (GWT), viz. 
TT = UmgNo To = {0} and 

T„ := {vi ...v„ e V|vi e T„_i and 1 < v„ < } 

(using the common tree notation vi...v„ for (vi,v„)). Put '■= #T„ for n € No- 
Further, let Zy denote the number of parasites in cell v € V and T* the set 
[number] of contaminated cells in generation n, i.e. 

T^={vGTn:Zy> 0 } and ( 1 ) 


for n e No- Over all generations, parasites sitting in different cells are assumed to 
multiply and share offspring into daughter cells in an iid manner. For parasites sitting 
in the same cell, however, the reproduction and sharing is conditionally iid when 
given the number of daughter cells. Formally, we let for each k gN 






/ e N, V e V, 


be iid copies of the Ng-valued random vector := and in¬ 

terpret as the number of offspring of the parasite in cell v which is shared 
into daugher cell vj given that cell v has k daughter cells. Assuming that initially 
there is one parasite sitting in 0, the numbers Zy of parasites in cell v are recursively 
determined by Z0 := 1 and 


^ 


= E = E^ffV G V,j G N, 


k>j 


( 2 ) 


where 


:= 0 is stipulated whenever j > k. Based on these variables, the pair 


(Tn, (•Zy)y(=Tn)n>0 


(3) 


is our branching within branching process (BwBP). The process of parasites is de¬ 
fined by 

:= E n G No, 

vgT„ 

and we further put 


7 := and pi k ■= for 1 <l < k. 


It will be a standing assumption throughout that the considered population starts 
with a single cell containing a single parasite. Nevertheless, it will sometimes be 
necessary to condition on a general number z G No of parasites in the ancestor cell. 
This will be expressed by writing P^, i.e. 
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Z^=z) = l, 

with corresponding expectation The index is omitted if z = 1, thus P = Pi and 

E = Ei. 

To rule out trivial cases, we also assume as in ||4j that 

0<7<oo, (Al) 

PI = P(A = 1) < 1 and P(^i = 1) < 1, (A2) 

7 ^ 1) > 0 for at least one l<j<k. (A3) 

By further assuming 

P2(^r>2)>0, (A4) 

we rule out the situation when all parasites in a cell share their offspring into one 
and the same daughter cell. In that latter case, the BwBP forms a branching process 
in iid random environment (BPRE) (see HI), for which the results to be derived here 
may be found in the literature, see e.g. Q m |2^ |22l |2^. By Thms. 3.1 and 3.2 
in |31, ( |A4| i ensures the extinction-explosion dichotomy for (^*)n>o and {2fn)n>o, 
i.e. these processes tend to infinity a.s. if parasites survive. Let Surv be the event 
of parasite survival and Ext its complement. Throughout this paper it is always 
assumed that parasites survive with positive probability, thus 

P(Surv) > 0. (A5) 

In combination with ( |A4| i and putting P* := P( |Surv), this implies v > 1 and 
t °o) = —7 oo) = 1, see ||4] Theorem 3.3]. 

Let us finally recall the definition of the associated branching process in random 
environment (ABPRE) (Z')„>o which forms a BPRE with Zq := 1 and an iid envi¬ 
ronmental sequence A := (A„)„>o taking values in {.if(A(-'^)|l < y < k < oo} (a 
countable set) with 

p(Ao = .if(A(>’^))) = ^ 

for all 1 < j < k < oo. Let (s) be the generating function of A„. The ABPRE de¬ 
scribes the evolution of parasites along a randomly chosen cell line (spine) through 
the tree, and Z' gives the number of parasites in the spinal cell in generation n, see 
in Section 2] for further details. This process is one of the major tools to study the 
BwBP, an important relation being 

P(z; >0) = (4) 

for all n G No, see 0 Proposition 2.2]. 

As pointed out in 0, the BwBP can be viewed as a multi type branching process 
with infinitely many types and comprises the process of type-A cells studied in 13 . 
Eor a more detailed account of the connections to other multitype branching models 
and related literature we again refer to 0. 
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2 Main Results 

In the following, we focus on the problem of finding the proper normalizations of 
the processes {■%*)n>Q and {3fn)n>o so as to obtain a.s. limits which are positive 
with positive probability. This problem has been studied in many branching models, 
see e.g. 0 El HI 112 HD m, and our goal is to provide analogous results for the 
BwBR We hrst concentrate on the process of contaminated cells (^*)«>o before 
turning to the process of parasites (.^)„>o. To prove the results for the latter, we 
use a spinal approach different from the one leading to the ABPRE and thus need 
more preparations and the construction of a size-biased parasitic branching within 
the branching cell tree. Let us stress once more that we only consider the case of 
a single ancestor cell hosting a single parasite, but that all subsequent results are 
easily generalized to arbitrary initial populations. 

Let A denote an arbitrary random variable with distribution {pk)k>o and (^n)„>o 
be the hltration, dehned by := { 0 ,f 2 } and 



for n > 1 


and '■= <y {\Jn>o ’^n) ■ It is obvious by dehnition that and are inde¬ 

pendent for all n > 0 , |v| > n and / > 1 . 


2.1 Growth rate of{^*)n>o 


Recall that {^n)n>o forms a standard GWR The Kesten-Stigum theorem El ensures 
that in the supercritical regime v" = is the right normalization of ,% in the 
sense that converges a.s. to a positive limit on Surv iff EN\ogN < How¬ 

ever, in order for this to be true also with 3^* instead of the parasite population 
evolving along a random cell line must have a positive chance to survive. In other 
words, the ABPRE must survive with positive probability. 

Theorem 2.1 {v " is a nonnegative supermartingale with respect to the 
filtration {^n)n>Q therefore almost surely convergent to an integrable random 
variable L as n ^ o°. Furthermore, 

(a) L — 0 a.s. iff one of the following conditions hold true: 

(i) V < 1. 

(ii) EN\ogN = °°. 

(Hi) Elogg^^(l) < 0 or Elog^(l-gAo(O)) =°°- 

(b) P(L = 0) < 1 implies {L = 0} = Ext a.s. 

The next two results address the question of growth rate in the case when L = Q 
a.s. Recalling v" = E.5^*/P(Z' > 0) from 0, the previous theorem tells us that 
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as n —?► oo should behave like its mean modulo an adjustment depending on 
the ABPRE. Since the environmental sequence of the ABPRE takes values in a 
countable set, EOl Theorem 1.1] states 

limP(z; >0)'/" = ,EA(1)“ =: P (5) 

where 

r 1, ifEiog4(i)>o, 

P = I v-'r, ifElog4(l)<0 and E4(l)log4(l) < 0, 

[ 1 A (' 7), otherwise. 


Hence, we may expect that the number of contaminated cells grows approximately 
like (vp)" and that a proper norming should not differ much from this sequence. 

Theorem 2.2 Suppose P(Surv) > 0, thus particularly v > 1. Then 
lim-log,5^* = logvp 

n^oo n 


If the ABPRE survives with positive probability, {^„*)„>o has nearly the same 
growth rate as the GWP {.%i)n>a (see Theorem 2.11, whence the Heyde-Seneta 
norming of {ST^)n>Q gives the right normalization for the process of contaminated 
cells in this case as well. 


Theorem 2.3 //'Elogg^(l) > 0 and Elog (1 —gAo(O)) < there exists a 

sequence (c„)„>o in (0,°°) such that c„+ilc„ —)■ V and c^^ £T* converges a.s. to a fi¬ 
nite random variable L* satisfying P(L* > 0) = P(Surv). Furthermore, the sequence 
{cn)n>0 is a proper Heyde-Seneta norming for (,%)n>0 as well. 


2.2 Growth rates of {2fn)n>o 

Recalling P(Surv) > 0, it is readily seen that E^ = 7" for all n > 0 and that the 
normalized number of parasites process 

IT„:=7-"ir„, n>0, 

forms a non-negative martingale with respect to h hence converges a.s. 

to an integrable random variable W. The following results show that (W„)„>o has 
very similar properties as a normalized supercritical GWP. On the other hand, it 
turns out that in order for W to be positive on Surv an additional condition besides 
E.^ilogiFi < 00 is needed, which guarantees that the partitioning of the parasite 
offspring into the daughter cells is sufficiently uniform. 
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Before stating the results, let us mention a related but weaker one by Biggins and 
Kyprianou M Prop. 8.1] on normalized multi-type branching processes in a very 
general setting comprising our BwBP. 

Theorem 2.4 The following statements are equivalent: 

(a) P(VP>0)>0. 

(b) EVP = 1. 

(c) {Wn)n>0 is uniformly integrable. 

(d) E(sup„>oW„) <oo. 

Theorem 2.5 The expectation ofW is either 0 or 1, and 


EW=l 


iff 


EJllog^l <oo 


and E 




< 0 . 


in which case P(VP > 0) = P(Surv). 


Our third theorem asserts that still grows like its expected value 7" on a 
logarithmic scale if log iFi = Thus, a proper normalization should not differ 
much from this sequence. 


Theorem 2.6 ^E 




< 0, then Wn 


1 jn 


1 a.s. on Surv as n 


The proofs of the stated result, especially Theorem 2.5 will make use of the 
size-biasing technique, which since the work by Lyons et al. m has become a 
standard technique in the study of branching models, see e.g. iKHEslElKIl 
mum ED and also 1^ for a similar construction in the context of multiplicative 
cascades. In Section]^ we will dehne a size-biased BwBP which is different from 
the ABPRE (Z^)„>o introduced in ID and in fact strongly related to a branching 
process in random environment with immigration (BPREI). The latter will therefore 
be discussed in the following section including the statement of limit results that 
will be useful for the analysis of the size-biased BwBP. 


3 The branching process in random environment with 
immigration 

The following can be seen as a stand-alone section of this article and does therefore 
not refer to the notation previously introduced. 

The Galton-Watson processes with immigration in hxed environment has been 
studied by many authors in the past, see a for the most important results and also 
references. In a multitype setting and random environment. Key ini and Roiter- 
shtein ll25l proved limit theorems in the subcritical case. Results for the single-type 
process in random environment for all three (subcritical, critical and supercritical) 
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regimes have been obtained more recently by Bansaye Eol. On the other hand, a 
theorem of Kesten-Stigum type for the BPREI, indispensable for our analysis of the 
BwBP, appears to be an open problem and is therefore presented below (Theorem 


HD- 


Turning to a model description, denote by the set of probability laws on No 
and by the subset of laws with finite mean. Let the environmental sequence 
= (‘^„)n>o = (^n,i) ‘^n. 2 )n >0 consist of fid random variables taking values in the 
set !^\ X endowed with the ff-field induced by the total variation metric. Given 
, let {Xn,k I (n, k) G No X N} and { |« G No } be conditionally independent families 

of iid No-valued random variables such that, for all n > 0 and k>\. 


€ •I'i') = and P(^„G-|‘i') = a.s. 


To ensure that immigration occurs with positive probability, we assume throughout 
this section that 


P(^o>0)>0. 


(6) 


The BPREI (Z„)„>o with environmental sequence ^ is then defined by Zo := 0 
and, recursively. 


z„ 

Z„+i ;= ^X„ k + (7) 

!=1 


for n >0. The Zjt „, k > I, provide the numbers of offspring of the individuals at 
generation n, while gives the number of immigrants at time n. It is clear by our 
assumptions that Z„ and {X^ > n,k > 1} are independent for all n > 0 which 

in turn ensures the Markov property for (Z„)„>o. Let 


the mean of As in the setting without immigration, we consider the super¬ 
critical case Elog/i.^g > 0, the critical case Elog/i^p = 0, and the subcritical case 
Elog/lf/o <0. 

Before stating the main results of this section, we recall the standard fact that 


Xn f 0: ifEXo<°o, 

limsup — = < (8) 

«->■“> ” loo, ifEX'o = °o. 

for any sequence (Z„)„>o of iid and non-negative random variables. 

The following martingle limit theorem of Kesten-Stigum type for the supercriti¬ 
cal BPREI will be of great use for our later analysis of the BwBP. The proof follows 
arguments of Asmussen and Hering in Q for the branching process with immigra¬ 
tion. 


Theorem 3.1 Let > 0 and recall that < oo a.s. 

(a) TfElog^ ^0 < oo, then there exists a finite random variable Z„ such that 
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lim 



— Zoo 


a.s. 


and the following assertions are equivalent: 

(i) F{Zoo > 0) = 1. 

(ii) F{Zoo > 0) > 0. 

(Hi) E((Xi,olog+Xi,o)/M'%) < °°- 

(b) ^Elog+^o = oo, then\ivasup^^„c^"Z„—ooa.s.foreverycG (0,°°). 
Proof (a) Defining the filtration 


ff (Zo,Zi,... ,Z„, {^k)k>o,‘^), « S No, 

thus ^0 = ’^{{^k)k>Q ,‘^)7 the sequence ((H^Jo P%V^^n)n>G is adapted and a.s. a 
L*-bounded, nonnegative and thus a.s. convergent submartingale with respect to the 
conditional measure P(-|^o) ^s the subsequent arguments show. We have 


z z 

]E(Z„+i|^„) = ^ E(Z„ + ^„ > ^E,{X„k\‘^) = 'Z^n a.s. 

k=\ k=\ 


and thereby 


E(Z„+i|^o) — E (E(Z„+i |^„)|^o) 

= e( ^E(Z„,,|^„) + ^„ 

\k=\ 


^0 


— IE(Z„|^o)/r^^+ ^„ a.s. 
for « > 0. It then follows by iteration that, for all n>Q, 


E 


Zn-\-l 


=^0 = E 


Z„ 

ok 


=^0 + ™ 






= y — < y — 

A;=0 n!=0 Pf/i k>0 n,=0 

< ^exp (log+^,t-El°gMf^ ] 


Ar>0 


;=0 


= L 

t :>0 \ 


k + 


k 

- (log+^,t-l^logM-2^; 
t ' !=0 


i:+l 


(9) 


a.s. (10) 


Since (^„)„>o and (/i^„)„>o are iid sequences and Elog^ < °°, ® and the strong 
law of large numbers provide us with 
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Elogn^p/g < 0 a.s. 


and thus the almost sure finiteness of the sums in (|^ and ( [TOl i. As a consequence, 
„>0 is indeed a -bounded and thus a.s. convergent submartin¬ 
gale under P(-|,^o) which leaves us with a proof of the equivalence of (i)-(iii). 
Denote by (Z„)„>o a BPRE starting with one ancestor, environmental sequence 
and no immigration. By ll2^ Theorem 2], (Z„/EZ„)„>o converges a.s. to a limit 
W as n —> oo, which is nondegenerate, i.e. q{‘^) := P(1T = 0|'^) < 1 with positive 
probability, iff (iii) holds true, i.e. E{{Xi Q\og^X i q)/< °o. So it remains to 
verify the implications 

P(Zo„>0)>0 ^ P(^('^)<1)>0 ^ P(Z„>0) = 1. (11) 

We show the first one by contraposition and assume that ) = 1 a.s. Note that 


n 4 

Z„ = Y.L^k,n-kii), 
k=0i=l 


( 12 ) 


where Zii „-k{i) denotes the number of individuals in the (n — generation of a 
BPRE started with the immigrant at time k and with reproduction laws given by 
‘^k,\T^k+iA ■ ■ ■ (see ifTTll and recall Zq = 0). Moreover, the (Z^„(/))„>o for k >0 
and ! > 1 are conditionally independent given and 


P((Zi,„(0)„>oe •IWi>« = u) = P((z„)„>oe-|^=u) (13) 

for P('^ e •)-almost all u € x As {‘^n)n>k = ^ and thus q{{‘^n)n>k) = 
q(^) = 1 a.s., it follows that 


'Zk,n—ki}) _ 1 ^k,n—k{0 


0 a.s. 



as n —^ oo for all i > 1 and A: > 0. By using these facts in ( [T^ , we now infer for each 
m e N that 



n m —1 ,, 

;=0 


as n —oo, where Z^ is a copy of Zoc and independent of j). To see 

this, one should observe that 


P(Tm,„ G -l^) = P(Z„-,„ G ■\{^k)k>m) 
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and the independence of and (^o, for any m = 1,The distribu¬ 

tional equation just derived, viz. 


Zoo 


d 


zL 



for all m G N, in combination with Zo„ < oo a.s. and 


m— 1 / m— 1 \ 

n = exp ^ log^^^ ^ °o a.s. 
j=0 \j=0 J 

by the strong law of large numbers obviously entails Zoo = 0 a.s. 

For the second implication in ([n) suppose now ) < 1) > 0, which partic¬ 

ularly implies Zn ^ with positive probability, see ll26ll . But this in combination 
with ([nil and (El easily implies Z„^ a.s. Fix £ > 0 and choose rj > 0 such 
that 

¥iq{^) <l-ri)> 1-e. 

For any A: S N, we then hnd that 


P(Zoo = 0|Z,,'^) = 


z„ 


= 0 


< 


lim , 


lim£ 

n—\oo I 1 


Zk,^ 


V .^in: 

= v{w=Qm)i>k) 


= 0 


k 


Zk.^ 


a.s. 


where Z„{i) describes the offspring in generation n stemming from the indi¬ 
vidual in generation k and thus behaves like the BPRE Z„ (modulo a k-shift of the 
environment). Since the population of the BPREI explodes almost surely and 
consists of iid random variables, we finally conclude 


P(Z„ = 0) < Eq{{^i)i>kf'‘ < E(1-TJ)Z*+£ £ 


which proves (a) because £ > 0 was arbitrarily chosen. 

(b) Let c ^ 0 3.nd notice tli3,t ^ 3,.s. for 3,11 n ^ Oby 0. Now use (j^ in 

combination with E log^ = oo to conclude 


lim sup — 

n^oo C 


> lim sup 

n^oo 




lim sup 




n 

= oo a.s. 


as claimed. □ 

As a consequence of the proof of the above theorem, we note the following corol¬ 
lary. 
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Corollary 3.2 In all three regimes, it is true that 

limsup —E (Z„|.^o) = 0 a.s. 

n—>oo C 


for each c > 1 such that Elog/i^.^^ < logc, where = c7((i*„)„>o,‘^). 
Proof. Let c > 1. By (|^, we have 


1 

c”+' 




< 




a.s. 


IfElog^^,,. >0, the proof of Theorem pJl has already shown that the sum on the 


right side is almost surely finite. Since the jXoK^, n £ No, are iid, we further get 


limsup]~[ = limsupexp j ^ log^^ ) ~ ® 

Ar=0 ^ \i:=0 ^ J 

by the law of large numbers, hence 

l™sup-^E(Z„+i|,^o) =0 a.s. 

n^oo C 

IfElog^,^^^ < 0, then the assertion follows by a simple stochastic comparison argu¬ 
ment (replace fZ by fZ' = ('^„')„>o satisfying Elog/r^j G (0,logc) and use that the 
assertion is true in the supercritical case). We omit further details. □ 


4 Size-biased branching within branching tree 

Unlike the size-biased construction used in HI to define the ABPRE, which was 
purely based on the cell tree, the following size-biased version BT of the whole 
BwBP will be obtained by picking a spine (random line) of parasites. Yet, since the 
spinal parasites are hosted by unique cells, this will again determine a random cell 
line as well (see Fig. [^l, but its properties are different from those of the cell line 
related to the ABPRE. 


Construction of the size-biased process 

Let 

n>0, 

be iid copies of the random vector independent of 

and (Av)veV^ where T,C take values in N and := is a 
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vector of random length T. These random variables have the the following distribu¬ 
tions: For A; S N, jc = (xi,... ,Xk) € Nq and 1 < m < we have 



P(f-k) - ^^^5=1 

(14) 


F{X^‘’^') =x\f = k) = =x), 

(15) 

and 

P(C = mX(*’^) =x,f = k) = , 

(16) 

i.e., C is uniformly distributed on {1,..given = x and T 

particular, 

=x,f = k,C = m) = yP(x(*’^) =x). 

= k. In 

(17) 


These random variables determine a random path (spine) through the parasites 
as depicted in Fig.[T]by the following 3-step procedure: 


Step 1 The root cell Vb = 0 splits into Tq daughter cells. 

Step 2 If 7b = k, = {xi , ...^Xk) and x = 'E!j=iXj, then the parasite in 0 has 
X descendants of which Xj go into daughter cell j for j = 

Step 3 The spinal parasite of the hrst generation is picked by Cq uniformly at 
random from the x offspring parasites and the cell hosting it is the spinal 
cell Vi of the hrst generation. 

The procedure is then successively applied to the spinal cell and its spinal parasite 
of generation n = 1,2,... So, being at generation n, the spinal cell Vb splits into 

Tn daughter cells, provides the offspring numbers of the associated spinal 

parasite and C„ the spinal parasite in generation n+\. All remaining parasites in I 4 
multiply independently with the law of We thus obtain a random cell line 

(V„)„>o with Vo := 0 and 

Fn+l := ViJJn 

for n > 0, where U„ denotes the daughter cell containing the spinal parasite of the 
next generation. Since the , n > 0, are iid, so are the Un, n > 0, and 

we get from 

P(fo =KUo = j) = for 1 < ;• < k < -. (18) 

7 

All parasites and cells not in the spine reproduce with the usual law. Thus, Z 0 = 1, 
and the children of each cell and their parasites in the size-biased BwBP 

, (F„, C„), (Zv) J 
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Fig. 1 A typical realization of a size-biased BwBP. Spinal parasites are shown as ■ and cells 
hosting these pai'asites as □. All other parasites and cells ai'e shown as • and 0> respectively. 
Bended edges are used to indicate the line of cells containing the spinal parasites. 


are given by 


- 

\fn, 

< 

II 


K = 




[n., 

ifv 


> 

1 

“r , 

II 

> 

Zv, •= j 

If: 

•J 

= 1^?,V 5 

ifvy^y„. 


for V G V with |v| = n and j G N. Finally, let T and ^ have the obvious meaning. 

It is important to point out that not only the spinal parasites but also those sharing 
a cell with them have a different offspring distribution as parasites sitting in regu¬ 
lar cells. This is caused by the fact that a spinal cell always produces at least one 
daughter cell wheras regular ones may die. The next lemma provides us with the 
reproduction distribution of a spinal cell and the parasites it contains. 

Lemma 4.1 The conditional distribution of (Tn,U„, (Zp ) given Zp , = 

the number of parasites in 14 - 1 , equals the distribution of {Tq^Uq^ under 

Pj/or all G N, and 










14 


Gerold Alsmeyer and Sdren Grottrup 


P;(7o = k,Uo^ I, Zj = Zjfor j = l,...,k) 


PkZl 

ZJ 


"^zi2j=Zjfor j =l,...,k). 


for allk&'H, 1 <l < k, and {zj)i<j<k & Nq. 

Proof. Let k gN, I G {1, {zf)\<j<k € Nq and z S N. Using the independence 

and i.A;> i.vgVj we then obtain 


i{Tq — k, Uq — f (Zj)i<y'<i — {zj)i<j<k) 


z-l 


To = k,Uo = I, = U- for j=h---.k 


i=l 


= L = (xi,... ,Xk), To =k,Uo = 1) 


Xi<Zi 


Z-l 


X P I = Zj-Xj for j=l,...,k 


Pk 


J=\ 

^z-l 


^ P I =Zj-Xj for j = l,...,k 


!=1 




rU,k) 


U,k) 


y Xj<Zj \!=1 


P^ HZ I vikP} 


= —EX 


7 


z,ei 




E^i.0 =2ifor7 = Pj(Z;=zyfor; = 1,...,^) 

Z 


=ZiWz{Zj=Zjfor j=\,...,k), 


= 1 


where ([T7J was used for the third equality. Since a random walk {S„)n>o with So = 0 
and iid increments Xi ,X 2 ,... satisfies E(Xi |5„) = S„/n a.s., we conclude the desired 
result. □ 


Dichotomy of the size-hiased process 

The next (common) step is to establish, by drawing on a measure-theoretic result due 
to Durrett M, equivalent conditions on the size-biased BwBP for the martingale 
limit W = lim„^„o Y ''^n to be finite. In the following, the dagger symbol t is used 
for formal convenience to indicate that a node v S V is absent in the cell tree and 
thus called a dead cell. Put ^ := Nq U {t} and define 

Zv := ZvljvGT} Ptl{v^T} ^rid Zy := Zvll|^gjj + 
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Then Zy = 0 means that v is a living cell with no parasites, whereas Zy = t indicates 
that V is not present in T and thus a dead cell. We put 

BT ■= (Zy)vgv and BT := (Zy)ygv 

and call these objects the branching within branching tree (BwBT) and the size- 
biased BwBT, respectively. It is important to note that previously introduced ran¬ 
dom variables of the BwBP and its size-biased counterpart, in particular ^n, can be 
defined as measurable functions of BT or BT. 

Let S := be the set of BwBT’s which assigns a nonnegative integer or t to 
each node of V. Put V<„ := {v : |v| < n} and V„ := {v e V : |v| = n} for n G Nq. Let 
be the standard (7-field on S generated by the projections on the components, and 
let C denote the sub-a-field which is induced by the projections on compo¬ 
nents in V<„. Clearly, {S^n)n>o is a filtration in and BT and BT are both (S, S^)- 
measurable. For n G Nq, we denote by BT„ and BTn the BwBT and size-biased 
BwBT up to level n, respectively, i.e. BT^ := (Zv)vev<„ ™cl BT„ := (Zy)ygv^^. 
Further, let Zn : § —> No denote the number of parasites in the generation of a 
host-parasite tree, thence = Zn{BT) and ^ = Zn{BT), and put 

w„ : S^> [0,oo), w„(t) := ^Zn(T) (19) 

for n G No. We further set w := limsup„^„w„. Thus w„ is ,5^„-measurable by defi¬ 
nition, and we have the representations 

Bn = w„ o BT and = w„ o BT. 

As a consequence of the following lemma, the uniform integrability of (W„)„>o is 
directly linked to the almost sure finiteness of W. 

Lemma 4.2 

(a) For all n G No, z = (zy)yGV<n € and u G V„ 

P(BT„=z,y„ = u) = ^P(BT„=z). 

(b) Let n G No and h : (S,,5^„) -G (M,^) be a measurable and non-negative (or 
bounded) function. Then, 

E(/z(BT)) = E{W„h{m)), 

in particular, for A G 

Q(A) = E(W„1{bxga}) = Qidt), 

where Qf) := P (BT G •) and Qf) := P(BT G •). 


(20) 
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(c) The following dichotomy holds true: 

(i) P(# < oo) = 1 EVK = 1, 

(ii) P(# = oo) = 1 P(VK = 0) = 1. 

Proof, (a) Since the statement for n = 0 follows from our definitions, let n G N and 
u = u'u„ for some u' G V„_i and u„ G N. Then by induction, the branching property 
we get for each (zv)vGV<n ^ 


and Lemma 


4.1 


P(BJ'„ — (Zv)vGV<ni^n — 

= P(Br„_l = (Zv)vGV<n_l,V„-l = u') 

X P((Zv)vGVn ~ (2v)vGVnJ^« “ = (Zv)vG¥<n-ll^n-l “ ^ ) 


— ~ (^v)ve¥<n-i) ri Pzv ((^v')v'gN — (^w')v'gn) 


jv|=«—I. vt^u' 




= ^P {BTn-l = (Zv)vg¥<._ J n ((ZvO = (Zvv')v'gn) 

t |v|=rt—l.VT^u' 


^ ((Zy') — (-^u'vOv'gn) 

^u'7 


= ^P(fi7’« = (^v)vG¥<J. 

(b) It suffices to show ( |20l i. Summation over u G V« in (a) yields 

^|u|=a 


’ {bT„ — (Zv)vG¥<n^ 


ytl 


' {BTn — (Zv)vG¥<n) 


for all (zv)vG¥<n ^ Therefore, we infer for all A G 


Q{A) = P(^GA) = P(^„GAn^^^") 

^|u|=n^u 


L 

L 


7” 

^|u|=«^u 


BT„ G t/(Zv)vG¥<n) 


yn 


P(Br G t/(zv)ve¥) 


— dQ — E(W„l{B7’g,4}). 
(c) Part (b) and fM] Theorem 5.3.3] imply for all A G 

2(A) = [ wdQ + Q{An{w = oo})^ 


which in turn leads to 
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ElV = [ w dQ = 1 — 2 (w = °o). 

Js 

The asserted dichotomy now follows. □ 

Remark 4.3 Since log+(wi(-)) is .5^1-^-measurable and nonnegative, the above 
theorem provides us with Elog+ Wi = EVTi log^ Wi, which in turn yields 

EJ^ilogJi<oo iff Elog.^i<oo. 


The process of parasites along the spine 


We take a closer look at the process (Zp )„>o of parasites along the spine and its 
recursive formula 


= 

Vn+\ 


E 

i=l 




(U„,Tn) 


-x;, 


(Un,fn) 


n e No- 


( 21 ) 


Observe that all but the spinal parasite in a spinal cell multiply with the same distri¬ 
bution, whereas the spinal parasite produces offspring according to a different law. 
We can figure the spinal parasite to be outside the cell and its progeny as immigrants 
of the spinal cell of the next generation. Then all remaining parasites in the spinal 
cell reproduce with the same distribution and we thus see that (Zp )„>o behaves like 
a BPREI. 

Theorem 4.4 Let (Z')„>o be a BPREI in iid random environment A = {A„)„>o tak¬ 
ing values in {{X^bk),k) _ = {j,k)) . 1 j }c ^ oo^ dvid such 

that 

- l)|(t7,f) = (;,k))) = 

for all \ < i <k <°° and Zq = Zp^. Then the law of — l),j>o equals the law of 
the BPREI 

Proof It suffices to verify that, for each n > 0, the conditional laws of Z'^j given 
and Zp j given Zp ,4,, coincide, for both sequences are Markovian. It 
follows from the definition of (Z')„>o and its recursive structure (see Q) that 


E \z’„=z,An=^{{X^'’’'^\x^^''^^ - mUJ) = {j,k)) 


= E 




fUP-i 


iU,T) = U,k) 


for all n,z S No and I < j < k < °o. Further, ( |2l] i and the iid assumption of the 
involved random variables implies 


A+r 


zp - 1 =z,{Un,T„) = ij,k) 


E 
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= E 


= E 


Ay(j] 






{UnX) = U,k) 

{0,f) = U,k 


for all 1 < j < A; < oo, hence 


E ( sX+i ' 


— 1 = z j = E 


Z=z 


for all n G No and z G Nq. □ 

We call the BPREI (Z')„>o in iid random environment A the associated branch¬ 
ing process in random environment with immigration (ABPREI) and denote by 

^/l„('S)= 52 ) ^{4„=if((A0A),f(t,*)-l)|(£/,f)=(j,i:))} 

j<k 

the generating function of the first marginal distribution given by A„. The process 
is called supercritical, critical or subcritical if Elogg^^(l) > 0,= 0 or < 0, respec¬ 
tively. 

Remark 4.5 There is a strong connection between the behaviors of the ABPRE and 
the ABPREI. Namely, if pj^ 1 and pk > 0 for at least one pair {j,k), i<j<k, 
then (see IH Section 2] for the definitions of different subcritical subregimes of the 
ABPRE) 

{ subcritical, | strongly subcritical, 

critical, iff ABPRE < intermediate subcritical, 

supercritical, [weakly subcritical or non-subcritical. 


This can be easily assessed by a look at the equation 

]Elogg'4^(l) = 52 (22) 

where we refer to H for the definition of the generating function gAoX- Since 
the function x i—.rlogx is strictly convex and ^ 1 with positive probability, 

Jensen’s inequality provides us with 

E4o(l)log4o(l) > %X(l)logE4o(l) > 

which combined with (|2^ shows the assertion. 
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5 Proofs of the main results 

Proof (of Theorem \2.1\ . Recalling the definition of (=^„)„>o from Section 1^ and 
noting the independence of and for each v S T„, the supermartingale property 
and thus a.s. convergence to an integrable random variable L of (v^"^*)„>o follows 
from 



< ^ E(iVv ^n) = L 


a.s. 



for each n > 0. 

If V > 1 and EA^logA^ < then (v ^^f)n>o is even uniformly integrable be¬ 
cause the obvious majorant £f,)n>o froms a normalized supercritical GWP sat¬ 
isfying the (Z log Z)-condition of the Kesten-Stigum theorem (see 0 Section 1.10]). 
Consequently, 



(23) 


where the second equality follows from 0. The theory of BPRE (see e.g. i) or E 
Prop. 2.3]) now implies in this case that L — Q a.s. if and only if condition (iii) holds 
true. If, on the other hand, V < 1, then < ^n=0 eventually, and if EA^ log = oo, 
then the Kesten-Stigum theorem implies 


L = lim —^ < lim — = 0 a.s. 

n —Voo n —^oo 


n^oo n^oo y‘ 


In both cases we obtain L = Q a.s., which completes the proof of (a), 
(b) Defining T„ = inf{m S N : > n}, we find that, for any n G N, 

P(L = 0) < P(L = 0 |t„ < °o)= °o) 



< P(v-'",^„:^0)” + P(t„ = oo) 


P(L = 0)" + P(t„=c«), 


+ E(t„ = °o) 
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where =5^ (v) denotes the number of contaminated cells in the generation of the 
subtree rooted in cell v € T*. Since P(L = 0) < 1 and Ext = {sup„>j ^n*< cxdJ- 3,.S. 
by Theorem 3.2(a) in our companion paper a, we arrive at the conclusion 

P(L = 0) < lim P(t„ = oo) = p I sup < oo ) = P(Ext), 

V«>i / 

which in combination with Ext C {L = 0} a.s. proves the assertion. □ 

Proof ( of Theorem \2.2\ . For each £ > 0, use Markov’s inequality it obtain 

^ \ V" 

I 


El 

n=0 


whence by the Borel-Cantelli lemma 




< 


lim sup 


( \ 

VE^„V 


1 /n 


< 1 a.s. 


But from Q, Q ^nd with Jensen’s inequality, we infer 

= vP(Z'>0)^/" ^ vp, 


as n —>• oo and thus ^ 

limsup-log,i5^* < logvp a.s. 

n^oo fl 

Left with the proof of liminf„^a, * log fPf > log vp a.s., we have —> oo a.s. 

on Surv by another appeal to 0 Theorem 3.2(a)]. Therefore, by Fatou’s lemma, 

oo = E lim inf 5;* < lim inf E ,5;*. (24) 

>oo n^oo 

We proceed with the following construction of a sequence (T*, „)n>o of sets of con¬ 
taminated cells for each m G N. 

Step 1 Put q := {0 } and suppose the root cell to host a single parasite. 

Step 2 Let j := T]], be the set of contaminated cells in generation m. From any 

of these cells, pick an arbitrary parasite and let j be the set of contam¬ 
inated cells in generation 2m containing at least one of their descendants. 

Step 3 Recursively define as in Step 2 with the help of T*, „ for each 

n = 2,3,... 

Then we obviously have 


:= #T;^„ P-a.s. 

for all n G No, and {Sm,n)n>o forms a simple GWP with offspring law P(=i5^ G •). It is 
supercritical for sufficiently large m by Let SurVm denote the set of survival of 
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{Sm,n)n>o, Obviously, Surv„, C Surv for all m G N. Fix mo such that P(SurVm(,) > 0 
and note that Surv^^ C Surv 2 mo ^ • C Surv a.s. because a GWP considered only at 

the points in time /Nq for any hxed 1 G N is also a GWP and survives if the original 
one does. Using these inclusions and the branching property of a GWP, we hnd 

P(Surv,^J = 

j>i 

> (i-P(Surv^j^)P(,;7,;^>5) 

for all s,kGN. Hence, 

P(Surv) > p( ySurVim^ ] = limP(Surv^^Q) > (1-P(Survy)^)P(Surv) 

V-t>o / 

for all s G N, and since P(Surv„,(,) > 0 is assumed, we arrive at 

y Surv^-mg = Surv a.s. (25) 

k>0 

by letting s tend to inhnity in the above inequality. 

Fixing now any m G moN and then k„ G N, /„ G {0,1} such that n — 
kntn + 1„, it follows on Surv^ that 

> E 


for all n G N, where ^ denotes the number of contaminated cells in generation n 
stemming from v G „. Using Jensen’s inequality, this yields on Surv^ 


= log+^* > -- Y. > log+5,„,i„ a.s., 

^m,kn veT* „ 


and the classical theory of GWP’s (for example the Heyde-Seneta theorem ||5] The¬ 
orem 5.1 in Chapter II]) provides us with 


lim inf - log > lim inf - log+ t 

n^oo Yi n^oo Yi ’ ” 

= —logE,^^ = logv -f —logP(Z^>0) a.s. 
m m 

on SurVm, where 0 has been used for the last equality. As m = kniQ for some k gN, 
we hnally obtain by using 0 and recalling ( |25l ) that 

lim inf-log ,3^* > logv + limlogP(ZL > 0) = logvp a.s. 
n k^oo kniQ “ 


on Surv. This proves the theorem. □ 
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Proof (of Theorem \2. 3| ). W.l.o.g. we may assume that EA^logA^ = for otherwise 
(v")n>o provides a suitable norming sequence by Theorem 2.1 For each a > 0 such 
that EA^ 1 L|^<qJ > 1 , we define 


co(a):=a and c„(a) := c„_i (a)E for n > 1, (26) 

which is obviously uniquely determined by a. Let (c„)„>o be a particular choice 
which is kept fixed herafter. Recall that (^„)n>o is a supercritical GWP with re¬ 
production law ^{N) and mean v. The classical theory of GWP’s (see e.g. jSj 
Chapter II]) tells us that any (c„(a))„>o provides a suitable Heyde-Seneta norming 
for {^n)n >0 with 


lim c«+iW _ - V c„(a)P(A^ > c„(a)) = 0. (27) 

Cn{a) 


Moreover, all (c„(a))„>o are asymptotically equivalent in the sense that there exist 
constants ri(a) G ( 0 ,<^) such that Cniffjjcn —>■ ?](«) as n — 

For n € No and v G V with |v| = n, put 

■^v(a) := A^vl{iv„<c 4 a)}, (28) 

and let T*{a), and be the obvious variables in a BwBP with under¬ 

lying cell process given by (A^v(a))veV- Classical theory asserts that the process 
(c^* (a) 5 ^j(a))„> 0 , is a -bounded martingale (see e.g. 19] Theorem 2 on p. 9]). As 
in the proof of Theorem |2.1| we derive for n G No 

/N^(a) 

= L E ^ 1{Z„,>0} I 

veT*(a) \ j=l 

Cfi [Cl) 

a.s. Hence, (c^*(fl),^*(a))„>o forms a positive supermartingale with E5^*(fl) < 
Cn{a)ja, and since the obvious majorant {cf^{a)^n{a))n>o is -bounded, there is 
an almost surely finite random variable L*{a) such that 

—» L*{a) a.s. and in L* (29) 

Cn[a) 

as n —oo. The rest of the proof splits into several parts. 

(1) Convergence oe SPf Ic„ 

With calculations as in the proof of ifT^ Prop. 1], we obtain 

for some n G No) 

= Y. = ^1 *,..., ^;(a) ^ ^;) 

n> 1 
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< L L = ^n-l = k, X{a) ^ X) 

n> 1 k>0 

< Y,^iN>Cn-i{a))Y,kV{^*_i{a)=k) 

n >1 k>Q 

< ^ P(A^ > c„_i(a))E^*_i(a) 
n>l 

< - ^ c„_i(a)P(A^ > c„_i(a)) —?► 0 as a —oo, 

“ n>l 

where the convergence follows from Hence, by ( |2^ , we infer for almost every 
O S f2 the existence of an uq such that for all a > aq 


Cfi 


Cnja) 3^*{a){(0) 

^nip) 


-4 ri{a)L*{a){ 0 }) 


(30) 


where rj{a) G (0,°°) should be recalled. Hence, c„ ^ converges a.s. to a random 
variable L*. 

(2) P(L* > 0) > 0 

In view of ( |30| , it suffices to verify P(L* (a) > 0) > 0 for some a > 0. Let A (a) = 
{A„{a))n>o be a sequence of independent random variables taking values in the set 
of probability measures on No such that 


F(A„{a)=^iX^J’^^)) 


Pk 

E(A^l{iv<c„{«)}) 


Cn{a) 

( \ Pk 

Cn+\\a) 


(31) 


for all n S No and Let further (Z^(a))„>o be a branching process 

in random environment A (a), and let gA„{a) denote the random generating function 
of the individuals in the generation. Recall that is the ABPRE with 

environmental sequence A. Clearly, 


P(z;(a) >0|A(a) =1) = P(z;>0|A=Z) 


P(Ao(a) = Ao,...,A„(fl) =1„) = —r^Y\P<k 

Cn\0) j .^0 

for any sequence of probability measures X = (A,),>o with Z, = AX ’and ji < 
ki < Ciia) for each i G No. Hence, by a straightforward adjustment of the calculations 
in the proof of Prop. 2.2 in Q, we obtain 

P(z;(a)> 0 ) = c^'(a)E,^*(a) 

for any n G No, and since c,7' (a) —>■ L*{a) in L*, we obtain 


EL* (a) = limP(Z'(fl) >0). 

n^oo 


( 32 ) 
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m —1 

8 Xn,i^) = E >m) 

k=Q 

be the generating function of the truncated random variable A m. As truncation 
reduces the reproduction, obviously 

EL* (a) = lim P(Z'(a) > 0) > lim P(Z'^(a) > 0), 

n—>oo n^oo ’ 

where (Z' (a))„>o is the branching process with environmental sequence A (a) and 

truncated reproduction laws. The truncation further ensures 

< °° a.s. 

n>0 

whence, by Theorem 1 in m, 

CO /n+l 

limP(z;^(fl)>0)>0 if E n 4 

n=0 \i=0 

Due to the given assumptions. Theorem 2.1 in ID provides us with the existence 
of a constant m > 0 such that 




( 1 ) 


< oo a.s. 


(33) 


0 < Eiogg;^^_,„(i) < 

A look at ( [3 T] i shows that 

p(A„(a)=.i^(zW))) = 

V / c„+](a) 

so that, by making use of (|27ll. 


n^oo pi^ 

, ^Pk -5- — 

c„+i(a) V 


P (Ao = if 


limElogg^;(„)^(l) = lim E pA:logE(z(-'’^^ A w) 

= E — logE(x(^’^(Am) 

\<j<k<oo ^ 

= ]Eiogg;^^,„(i). 


Moreover, for any x >0, 

P(log*4„(«),m(l)>^) 


Cn{a) av 

kPk < 


Pk 


l<;<t:<c„(fl).l<;<A:<co, 


log 


av 

Cl (a) 


log*M;,/t>^ 


P(log*gX,m(l) >-1^) 







Blanching within branching II: Limit theorems 


25 


and therefore an extension of the law of large numbers (see im Theorem 2.19]) 
ensures the existence of an a.s. hnite random variable M such that 


^ > ^Elogg;i^ ,„(l) > 0 for all n>M. 


k=o 


But from this, we hnally deduce 


«+i 1 

EBv 


M-ln+l 1 

En::r 


n >0 i =0 n =0 /=0 


n-\-\ ^ 

^exp( -^log4.(^),„(l) 


n=G 


M—\n+\ 1 

S En:7 


fi=0 i=0 n>M 


r=0 


exp --Elogg;\^,„(l) 


~\ «+l 


< oo a.s. 


and thereupon EL* (a) > 0 by an appeal to 


L* VANISHES ONLY ON Ext 


We adapt the proof of Theorem 2.1 b) and set T„ := inf{m S No|^ > n} for each 
n S N. Then 


P(L* = 0) < P(L* = 0|t„ < °o) +P(t„ = oo) 
lim AuL. ^ 

< Fi^:/c,n^0T+F{Tn=oo) 

= P(L*=0)" + P(t„=c«), 


< °o + P(T„ = oo) 


where entered in the penultimate inequality. As P(L* =0) < 1, the proof is 
completed by letting n tend to inhnity and an appeal to im Theorem 3.3]. □ 

Proof (ofTheorem \2.4\ . The implications “(iv) ^ (Hi) ^ (ii) ^ (i)” follow directly 
from standard martingale theory so that it remains to argue that P(VT > 0) > 0 im¬ 
plies E (sup„>o Wn) <°°. Modulo minor modihcations the subsequent proof follows 
the arguments of HB Lemma 2] and E Lemma 2.6 in Chapter II], and we estimate 
the tail probabilities of sup„>Q W,,. 

Let P(IT > 0) > 0. Assuming the existence of constants a,b > 0 such that 


P (W > at) > bF 


( sup W„ > t 

\n>0 


(34) 


for all t G it follows by a standard computation that 
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which proves the implication “(/) => (iv)”. So we must only verify p4| ) to complete 
the proof of the theorem. 

Proof of ([3^: Clearly, P(VP > 0) > 0 implies EVP > 0, and for each c G (0, EVP) 
we can find some w > c such that E(VP A w) > c. Fix f G [1, °°) and define 

En '■= \ ^n> t, sup Wk < t 
for n G No- Then, for any a > 0, 

P(VP>af) > P(vP>aA supVP„>f) = ^ P(VP > af|£„)P(£„). (35) 

V n>0 / „>0 

For V S V and « G Nq let denote the number of parasites in the n'* generation 
of the subtree with root cell v, the latter containing Zy parasites. Since (VP„)„>o is 
a martingale (under any P^), we obtain the almost sure convergence of ^ 

conditioned under Zy and denote its limit by VP^''^. For all n G No, we then have the 
representation 


1 




(v) 


1 


w = —lim y a.s. 

7" 7^ 7" yS* 


Consequently, 


P(VP > = P ( 4 E >at En] 

V veT* / 


^ y wi-) > — 


En 


> 


1 


E W^''> > a 


vgt: 


En 


> P E (W^^''^AZyW)>, 

V vgT* 


En 


¥{En)-^ / E (lV^''^AZyw)>a .^n] 

je„ y =z„ j 


dF. (36) 


For Z 0 = z G No let j denote the number of parasites in generation A: G No stem¬ 
ming from the ancestor parasite 7 G {1,..., z}. If any of these ancestors has at most 
wy^ in generation k, then the total number of offspring is at most zwy^, i.e. 
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7=1 

As a consequence. 




-a.s. 


V/=i 


1 


Azw) = Ej l^lim ^ j Azw 


1 


- ^H I™ ^■^>7 11 ^) ) = Aw) > zc 


for all z G No, which in turn implies 

El ^ L (^^''’aZ.w 


vgt; 


1 


J '^"vgT,*z=1 

^ i E = c a.s. 

vGT! 


for all n G Nq. Let us put W„{w) := y Lvgt; AZvw) and note that W'„(w) < w 
a.s. for all n G Nq. Then it follows for all a G (0,c) that 

rw 

c < E(W'„(w)|,F„) = / P(W'„(w) >x|,^„) dx 
Jo 

rw 

<a+ ¥ {Wn{w) > x\-^n) dx < a + (w — a)¥ {Wn{w) > a\,^n) 

J a 


and thus 


V{Wn{w) > a\,^n) > -—— a.s. 


w — a 

Plugging this inequality into ( [36l l with a := cjl and setting h := cji^w — c), we 
obtain 

P(lT>flf|£„) > b 

for all « G No and f G and in combination with (|T5]l this finally yields 


P(1T > flf) > ^ P(1T > flf |£„)P(£„) >hY^ P(£„) = /tP ( sup W'„ > f 


n>0 


n>Q 


n>0 


for all t G that is (|T4|i. □ 

Proof (of Theorem \2.5\ . First note that EVT G {0,1} is already verified by Theorem 
|2.4| The remaining proof is quite long and therefore split into several parts which 
are proved independently. 

Assertion 1; P(1T > 0) > 0 =a P(VT > 0) = P(Surv). 
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Proof of Assertion 1: Let T„ = inf{m G > n}. A similar argument as in the 

proof of Theorem |2.1[ h) yields 

P(1L = 0) < P(1L = 0)" + P(t„ = °o) 

for all n G N. Since P(VL = 0) < 1 and by ||4] Theorem 3.2], this further implies 
P(1T = 0) < limP(T„ = °o) = P ( sup < oo ) = P(Ext). 

Vm >0 / 

and then the assertion, for Ext C {W = 0}. □ 

Assertion 2: E^I log < oo and E (^^^log^^^ <0 E1T = 1. 

Proof of Assertion 2: To prove the stated result, we use the size-biased tree intro¬ 
duced in Section]^ and show that VT := limsup^^^W], is almost surely hnite. Then 
EIT = 1 follows by the dichotomy in Lemma [43] ]c). 

Recalling the notation of the size-biased process, we have the recursive represen¬ 
tation 


Av Zv 


^+1 = L = !:%,.+ I «>o. 

;=i vGf„\{v4-i=''=i 


Dehne the CJ-held 


Using the above recursive formula of ffn+i, we obtain 


Tn 


E{X+i\^) = E 


u=i 

= lE(E2i 


V„i 




+ E| E EE^if-’ 

^vGf„\{v,.}r=i'=i 


'VnJ 


Zv / Wv 


+ El E EE(Ex';'"’ 

\vGf„\{v„}'=i \r=i 


=7 




-f 


7e(^„|^) 


< ... < Er-^'E 


A:=0 




'Vki 


a.s. 


(37) 
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for all n >0. Using the definition of the size-biased variables and the fact that, for 
fixed {j,k), the / G N, v S V, are iid, we further obtain 


k=0 7=1 


xU,Tk) _,_]g 




V 


■ihTk) 


= x;r"-‘i;(sF’ 

k=0 7=1 


E ^ - 

1=1 


a.s. (38) 




Thus, letting n tend to inhnity on the right hand side, leads to 


E 




£ e4exF‘’ + e>(% 

^>0 ^ j=\ k>0 ' 



=-Ji 


='-h 


(39) 


a.s. for all n G Nq. Next is to show that both sums J\ and J 2 are almost surely hnite. 
Recall that 7 > 1 as a consequence of P(Ext) < 1 and ^ Theorem 3.3]. 

Finiteness oe Ji: By dehnition, the A: > 0, are iid with the same 

law as ^\. Moreover, ES]log.S°i < 00 is equivalent to ElogiFi < 0 ° (see Remark 
|4.3| l so that, by 


1 

lim - log 

it-ioo k 



= 0 a.s. 


which in turn entails that with probability one 


eventually, 

7=1 

in particular Ji < a.s. 

Finiteness oe J 2 : By Theorem 

vironment {U„,T„)„>o with immigration sequence {X„ — l)„>o- Consequently, 

^ , i G No, is the (random) reproduction mean of parasites in cell V, , and thus of 
the hrst marginal distribution of individuals in the generation of the ABPREI. As 
previously pointed out, E.^i log.S°i < 0 ° implies ElogS] < 0 °, and thus the immi¬ 
gration components satisfy 

Elog+(X(|^°’^“^-l) < Elog^ < 00. 


4.4 


(Zp — l)„>o is a BPREI in iid random en- 
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By the assumptions in the theorem and Eg^^(l) = 7 /v (see Q), we get 


E j = E4o(l)log4o(l)-^log7 < 0, 


whence, by an appeal to 


El°g^c/o,fo = < '°g^- 


(40) 


Consequently, Elog/r^^ = logc for some c € (1,7) and Corollary 3.2 applied to 
(Zp — l)„>o, provides us with 


Zoo := sup—E("zp — 1 | 


< °° a.s. 


We are thus led to a new upper bound for J 2 , namely 

k % 


-J Em, 


A:>0 


t=l 


■^j-Tk 


< Zoo ^ exp 


k>0 


C 1 


log- + -log+ I 




Tk 


a.s. 


Using Jensen’s inequality and ([T4li, we obtain 


7b 


Elog+ = EE(7b = ^)log+E 


(j.k) 


^7=1 


Ar>l 


V7=l 




k>\ 


W =1 / 

k 




A:>1 


v/=i 


7=1 


= -EiTilogiTi < 00, 
7 


and since the LJLi M ■ f , ^ > 0 , are iid, ([^ yields 


(41) 


limsup-log+ £M,fo = 0 a-s- 


, 7=1 


Hence, J 2 <°° a.s. follows when using this fact in (|4T]i. 
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Having verified that J\ and J 2 are a.s. finite, inequality ( [39] l gives 
supE(W„|i^) < 00 a.s. 

n>0 

while Fatou’s lemma ensures almost sure hniteness of liminf„^.,o i.e. 

P(liminfW„ < 00 ) = 2 (liminfw„ < 00 ) = 1 , 

where Q = V{BT £ •) and = w„ oBT. It remains to show that (w„)„>o converges 
2-a.s., because then W = liminf„^,x,VF„ and W is almost surely hnite, which com¬ 
pletes the proof of the theorem. 

We show that (l/w„)„>o is a Q-supermartingale with respect to the hltration 
(c5^„)„>o dehned in Section|^to which is adapted by dehnition. For each n > 0, we 
have 

g(w„ =0) = f w„dQ = 0 
J{Wn=0} 

by Lemma |4~2| h). Let E/> denote expectation with respect to a probability measure 

we infer for any A G Q =^«+i 


[ dQ = 

Ja ^ V>v«+i / 


< 


where the last equality follows by making the previous calculations backwards. We 
have thus verihed that (l/w„)„>o is indeed a Q-supermartingale with Eg(l/w„) < 

Eg(l/wo) = 1. It hence converges Q-a.s. by the martingale convergence theorem 
which completes the proof of the assertion. □ 

Assertion 3: Epilog^ 1 = 00 or e(^^^ log>0 ^ P(W = 0) = 1 
Proof of Assertion 3: First note, that our basic assumptions ( |A4| i and ( |A5| l imply 

P(^(i)e{7,0})<i, 

for otherwise jJLj k G { 0 , 7 } for all 1 < j < A: < 00 with pk> 0 and so 


1 


-dQ = E; 


JA 

Ee(— 

Vlt'n+l 
2(An{w„+i > 0}) 
Q{An{w„ > 0}) 

/ — dQ, 

Ja Wn 


Wn+l 


l{An{w„+i>0}} 


'^«+il{An{w„+i>0}} 


P. By using Lemma |4^ and Remark 4.3 


(42) 
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k 

7 = L L Pij.k = rL Pi#{l <j<k-. > 0) > 0}. 

k>a i=\ k>o 

'-V-' 

=;c = l 

Since c denotes the mean number of cells that are able to host parasites, we infer 
Ez(^i*) < c = 1 for all z G N. Hence, 


= E ^ < E,!?;; < ... < 1 . 

\veT*z=l / 


But this contradicts ii Theorem 3.3] and so ( |42l i must hold. 

Using again the size-biased tree defi ned i n Sectionj^ we show that P(1T = oo) = 1 
for IT := limsup^_ ^ 

First, note that 


W],. Then Lemma 4.2 c) gives P(1T = 0) = 1. 


Wn 


1^1 ^ 

y Zy > Tzg > 

1/71 yll 

vGf„ ‘ 


1 

yfl 


Tii-l 

;=i 


a.s. 


for n G N. Since E log 3fi=°° gives E log jFi = oo by Remark |43] and the random 

sums n G N, are independent and identically distributed as .Si, we 

infer 


^ 1 T ) 

limsupW„ > limsup— ^ 

n^»o 7" 


= limsup exp | -log V logy] = oo a.s. 

y i“i ^ ) 

by another appeal to (|^. This proves the assertion in the case when ESI log S] = < 
Suppose now ESI log S] < oo. Once again, by the definition of 1T„, we get 


1 V- 1 1 .'- 

W„ = — TZy > —Zp > —(Zp -1) a.s. 

yti i-j “ — yn V„ — yn'- Vi, > 

vgT„ 


(43) 


for n G Nq. As stated earlier, (Zp — l)„>o forms a BPREI in iid random environ- 

ment {Un,T„)„>o and with immigration sequence (Z„ — l)„>o- The assumption 

E(4o(l) G {7,0}) < 1 implies /Tyi ^ 7 for some ( 7 ,^:), j < k, with pk > 0 and 
P(Z(7U > 0 ) > 0 , thus 

P(Mc/o.fo ^ 7) > 0. 

Moreover, ES] log S) < 00 implies 


(44) 





Blanching within branching II: Limit theorems 


33 


Elog+- l) < 


as pointed out before. By adapting the argument in (|40|l, we find 


> logr > 0. 


(45) 


Hence, by Theorem |3.1 [ a), there exists an almost surely finite random variable Zc 
such that 



lim -- = Zoo a.s. 


(46) 


Theorem|3.1[a) further ensures that Zoo is a.s. positive, because 





^ l<j<k<oo 



where © has been utilized for the second equation. From ( |43| l, ( |44| ), ( |45] l, ( |46l l and 
the fact that the up f, t G No, are iid, we hnally obtain 



W = lim sup W„ > Zoo lim sup 



by the law of large numbers or, equivalently (Lemma [4.2[ c)), W = 0 a.s. □ 

This completes the proof of Theorem |2.5| □ 

Proof (of Theorem |2.6p . Since Wn converges a.s. to a finite random variable, it fol- 

\/n 

lows immediately that lim sup„_j,oo VTn < 1 a.s. 

To derive the lower bound, we distinguish two cases and use the truncation argu¬ 
ment given in ifl^ . Recall from (|42li in the previous proof that (|A4|i and (|A5|l imply 


^(4o(i)e{7,o})<i. 


Case I; Let N be bounded, i.e. N < c a.s. for some c G (0,°°). For a > 0 and 
I < j < k < c, we dehne 
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and let (^„(a))„>o denote the associated process of parasites and {^* {a))n>o the 
process of contaminated cells having the truncated reproductions laws. Put fur¬ 
ther 7 (a) := EjFi(a) and let gAo.ai^) and W„{a) have the obvious meaning. Since 
^„*ia) t as a —>■ for each n € N, we see that ¥ 2 (^ 1 *{a) > 2) > 0 as well as 
sup„>oE^*(a) > 1 for sufficiently large a > 0, which entails P(i?),(a) —)• 0) < 1 
for such a by 0 Theorem 3.3]. Moreover, 

E4„(l)log^ 

= E 4 ^_„(l)log 4 g_^(l) - Eg(^_^(l)log 7 (a) 

< Eg(^(l)log4o(l) - 

\ asa^oo 

= EgX(l)log^^ < 0, 

for 7 (a) increases in a. Therefore we can find ap > 0 such that 

< »■ 

for all a > ap. Since E^i(a)logiFi(a) < aclogac, Theorem |2.5| implies the exis¬ 
tence of a finite random variable VT(a) such that W„{a) —>■ W (a) in L* as n —>■ 00 . In 
particular, P(lT(a) > 0 ) > 0 . 

Now fix any e > 0 and a > ap such that 7 (a) > (1 — e) 7 . Then 
Ej;(a) = 7 (a)"> (l-e)" 7 " 

for all n > 0. Let {3f„ i;{a))n>o be the number of parasites process, when parasites 
produce offspring according to the original reproduction laws up to generation k 
and with the truncated laws from generation k+1 onwards. By the previously es¬ 
tablished lower bound of the means, this yields 

Eir„,,(a) = Y^^X-k{a) > (l-e)" 7 " 

for all k, n > 0 with k < n. Moreover, we find that 

> Xda) ^ ^n,k{a) > 1 y 
(l-e)" 7 « - (l-e)" 7 « " E^„,i(a) " 

where .^^(^[(a), v € T*, are iid random variables with the same law as ^_<-(a) 
when starting with a single parasite. Because of our choice of a, taking the limit in 
the above inequality yields 
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liminf 




n^oo (1 _ 


> liminf —r > ———^ 


= 4 E 


v€Tk 


where (a), v € T*, are independent and distributed as W(a). Recalling that 
is the (7-algebra of the A:-past, we get from this inequality 


P (liminf->0 >P( 




vgt? 


= l-V{W{a)=0y^k a.s. 

Now use P(iy (a) > 0) > 0 to obtain upon letting k tend to infinity 



Surv = {•SC^°°} C /liminf 7 —— -> o| 

^ " I (l-e)"7" J 

and then finally 

liminfWn'^" > 1 — e a.s. 

n^oo 

on the survival set Surv. This proves the theorem in the first case. 

Case II: Let N be unbounded. We use truncation to reduce to bounded N and 
make use of the results just shown for that case. For b>Q, put 


Let {2Yn{b))n>o be the associated number of parasites process and (=3^*(^))n>o the 
process of contaminated cells having the truncated reproductions law for the cells. 
Further, let Y{b) = EiFi (b), v{b) = (b) and gAo,b be the generating function of 

the ABPRE of the truncated BwBR For the truncated process, we have 

l<j<k<b 


as well as 

E4o,t>(l)log7(^) = = Eg;^^(l)log7 € (0,°o). 

Putting these two equations together and using yib) t 7 as — > oo, we obtain 

“ %X(l)log7 

= EgX(l)log^^ < 0. 
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Hence, for each e > 0, we can fix as in Case I some > 0 such that 

7 (fi)>(l-e) 7 , P2(^i*(fi)>2)>0, P(ir„(fi)-^ 0)<1 

and fo(l)log ^ ^ other words, all conditions of Case I are fulfilled, 
and so 


liminfwi^” > (1—e)liminf 


7(fi)« ) 


> 1-e 


a.s. 


This completes the proof. □ 


Glossary 


iPk)k>0 

T 

T„ 

N, 

% 

Zv 

BT 


i^j.k 

7 

V 

Wn 

{^n)n>Q 


N, 

Zv 

f„ 


Wn 

(Zn)n>0 


offspring distribution of the cell population 
cell tree in Ulam-Harris labeling 

subpopulation at time (generation) n [= {v € T : |v| = n}] 

number of daughter cells of cell v € T 

cell population size at time n [= #T„ = ^vgV -^v] 

number of parasites in cell v 

branching withing branching tree [= (Zv)vgv] 

number of parasites at time n [= ^vgT Zv] 

population of contaminated cells at time n [{v e Tp ; Zv > 0}] 

number of contaminated cells at time n [= #T*] 

given that the cell v has A: daughter cells v 1,..., vA:, the j component 

^/v of this Ng-valued random vector gives the number of offspring 

of the parasite in v which is shared into daughter cell v j. 

generic copy of the / G N, v € V with components X^^k) 

= EX^j’'^'> 

mean number of offspring per parasite [= E^i = Y.k> \ PkY!]=i Pj,k] 
mean number of daughter cells per cell [= EN ^LkyiPk] 
normalized number of parasites at time n [= Y oBT] 

infinite random cell line in V starting at Vq = 0, where Vn denotes 
the cell containing the spinal parasite picked at time n. 
number of daughter cells of cell v in the size-biased tree T 
number of parasites in cell v in the size-biased tree T 
number of daughter cells of the spinal cell Vn 
size-biased branching withing branching tree [= (Zv)vgv] 
number of parasites at time n under size-biasing [= Lvgf Zv] 

normalized number of parasites at time n [= oBT] 

associated branching process in iid random environment A — {A„)n>o 
with immigration (ABPREI) and a copy of (Zp — l)„>o- 
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